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ABSTRACT 

Context. Observations suggest that a relationship exists between the driving mechanism of roAp star pulsations and the heavy element 
distribution in these stars. 

Aims. We attempt to study the effects of local and global metallicity variations on the excitation mechanism of high order p-modes in 
A star models. 

/Wefhods. We developed stellar evolutionary models to describe magnetic A stars with different global metallicity or local metal 
accumulation profiles. These models were computed with CLES ("Code Liegeois d'evolution stellaire"), and the stability of our 
models was assessed with the non-adiabatic oscillation code MAD. 

Results. Our models reproduce the blue edge of the roAp star instability strip, but generate a red edge hotter than the observed one, 
regardless of metallicity. Surprisingly, we find that an increase in opacity inside the driving region can produce a lower amount of 
driving, which we refer to as the "inverse /c-mechanism". 
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1. Introduction : general context 

1.1. Tiie magnetic Ap stars 

The magnetic Ap stars are slowly rotating, chemically peculiar 
A and sometimes F type stars. They exhibit non-uniform dis- 
tributions of chemical elements, both laterally across their sur- 
faces and vertically with height in their atmospheres (see Sect. 
I2.2I 1. They have strong organised magnetic fields that are pre- 
dominantly dipolar, although evidence of deviations from dipo- 
lar structure are observed. The magnitudes of the fields range 
typically from a few hundred to a few thousand Gaus s (e.g . 
Landstreetl|T992l: iMathvs etaHUQQ?!: iHubrig etal.ll200l 120051: 



Wade et al. I2000h. They can h owever be str onger as d emon- 
strated by Hubrig et all ( |2005|) and Rvabchik ova et all ( |2006|) 
who detected magnetic fields stronger than 20kG in cool Ap 
stars. 



1.2. The roAp stars 

The coolest subgroup of magnetic Ap stars is the Ap SrCrEu 
group (^ 6400K-10000K). While most of these cool Ap stars do 
not show any oscillations, some exhibit single- or multi- period 
pulsations ranging from 4 to 21 min. These pulsating stars are 
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called "rapidly oscillating Ap stars" ( roAp stars) as opposed to 
non-oscillating stars called noAp stars. 

Both groups do not differ significantly in rotation and mag- 
netic field strength. The physical processes that create either 
roAp or noAp stars are not yet understood, although observers 
attempt to identify observational differences between roAp and 
noAp to solve this puzzle. We discuss this more extensively in 
Sect. 2. 

Until now 35 roAp stars have been detected (Q Curtz 198^ 
'Martinez ^aTII 199 1': 'Kurtz et alJll994tlMartinez &~Kurtz, 19941: 
Girish et al. [ I2OOU lElkin et al.l l2005h . Their rapid pulsations are 
interpreted as high-overtone, low degree, non-radial p-modes. 



2. Observational context 

2. 1 . Position in the HR diagram 

The roAp stars have effective temperatures of between 6400K 
and 8400K. They occupy mostly the main sequence part of 
the classical instability strip, in a similar way to 6 Scuti stars. 
However, a few are cooler than 6 Scuti stars: e.g. HP 213637 one 
of the coolest confirmed roAp stars (Tefj - 6400K, iKochukhovl 

Comparing obse rvational points to the models of 
Schaller et al. ( 19921) (wit h the solar in itial composi t ion of 
Anders & Grevessd (Il989l) and Z=0.02), iHubrig et all (l2000h 
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concluded that roAp stars are located on average slightly above 
the zero-age main-sequence, and significantly below the TAMS. 
This conclusion needs to be reconsidered however following 
the discovery of tw o apparently evol ved roAp stars: H P 213637 
(lKochukhovl2003l) and HD 116114 (.Elkin et al.ll2005l) (see Sect. 
l2.1.1l l.We demonstrate however in the following sections how 
the evolutionary stage of roAp stars depend strongly on their 
initial chemical composition and, in particular on their global 
metallicity. 



2.1.1. roAp versus noAp 

Figure [T| displays the effective temperature and the luminos- 
ity of roAp stars, noAp stars, and roAp candidates. The left 
panel shows the HR diagram of stars for which both the lu- 
minosity and the effective temperature are available in the lit- 
erature (19 roAp, 18 noAp, an d 4 roAp candidates tak en from 
iKochukhov & Bagnulol ( l2006l) and iNorth et alj ( Il997h and the 
right panel shows the effective temperature of magnetic Ap stars 
for which only the effective temperature is known (23 roAp, 4 
roAp candidates, and 24 noAp). 

From the left and right panels of Fig.[T] it is evident that roAp 
and noAp stars have a small overlap in effective temperature, the 
group of roAp stars being cooler than noAp stars. We also note 
the absence of cool non-pulsating Ap stars with effective tem- 
peratures lower than 7500 K, which suggests that all magnetic 
Ap stars up to a transition temperature of about 7500K-8100K 
may be affected by an excitation mechanism that is sufficiently 
strong to stimulate pulsations. 

Observations of the external parameters of roAp and noAp 
stars have suggested that noAp stars are m ore evolved and more 
luminous than roAp stars dNorth et all 1997 ; Hubrig et al. 2000). 
However, as discussed by 'Hubr ig et al.l (l2000i) : iCunhai (i2002l) : 
[Kurtz et al. (2004), this pattern could be due to observational bi- 
ases. According to iCunhal(l2002b . this bias could be related to the 
frequency of the excited oscillations, which are expected to be 
lower in more evolved and more luminous stars a nd possibly un- 
detectable. As mentioned bv lHubrig et al.l ( l2000l) and lKurtz et al.l 
(2004), most of the confirmed roAp stars were discove red by 
Martinez and Kurtz using photometric observations (e.g.lKurtz 
19821; iMartinez et al.|[T99lt iKurtz et allll994l;lMartinez & Kurtz 



19941) . In their studies, they preferentially searched for stars with 
negative Stromgren 6ci indices, because they are characteristic 
of strong chemical peculiarities in cool Ap stars and therefore 
provide good roAp star candidates. However, 6ci increases with 
luminosity so that a peculiar evolved cool Ap star may show a 
normal Sci. By focusing their surveys on stars with negative Sci, 
Martinez and Kurtz limited their research to main-sequence (or 
close to main-sequence) s tars, which i ntrodu ced a selection ef- 
fect into their discoveries. iKurtz et al.l (|2004|) also remarked that 
the highspeed photometric searches of Martinez and Kurtz were 
inadequate for determining periods of up to 20-25 min, and may 
instead introduce a bias into the determination of the roAp fre- 
quency range upper limit. 

Observations have indeed detected luminous and/or evolved 
roAp stars. From a comparison between the atmospheric pa- 
rameters of the cool roAp star HD 213637 and evolut ionary 
ti-acks computed by [Schallef etTali (1 1992b . iKochukhov! (|2003|) 
concluded that it is either located close to the end of its 
main-sequence lif e or has already ascended the giant branch. 
lElkin et aTl ( l2005b confirmed the existence of the luminous, 
evolved, lower frequency roAp stars, with the discovery of a 
21 minute pulsation in the cool magnetic Ap star HD 116114 
(Teff - 7413K, log(L/Lo) = 1.32). One implication of this dis- 



covery is that there may be more luminous roAp stars with yet 
longer periods. A search is now being undertaken to discover 
more of these stars and study them in detail. If they are found in 
significant numbers, then the luminosity difference between the 
roAp and the noAp stars may disappear, leaving an unanswered 
question: why do some Ap stars pulsate and do not others ? 

2.2. Chemical composition 

2.2.1. Surface Composition 

The cool magnetic Ap stars are chemically peculiar stars, which 
exhibit strongly anomalous, non-uniform element distributions. 
Their chemical peculiarities are distributed in patches about their 
surfaces, which are generally measured to contain underabun- 
dances in some light elements, such as C, O and sometimes N; 
almost solar to significantly low abundances of Fe and Ni (as 
well as other iron-peak element abundances); overabundances 
of Cr an d Co; and significant overabundances of rare earth 
elements (Kupka et al]|l996l; iGelb mann et al.lll997^ 'Gelbmannl 
1998; Ryabchikov a et al.lll997all b. 2000, 2005; Gelbmann et al] 
2000). An iron underabundance appears to be a common fea- 
ture of the chemical composition of the atmospheres of roAp 
stars with ef f ective temperatures of b elow 7000 K ( HD 21 3637, 
IKochukho v* ("2003); HD 101065, iCowlev et all JIOOOl) ; HD 
217522, lOelbmann (1998)), while, for hotter stars, iron appears 
to be close to the solar abundance values (IKochukhov! 120031 ; 
iRvabchikova et al.ll2004 . 

A Doppler imaging inversion technique has been developed 
(Piskunov & Rice' 1993) and used (IKochukhov et al.l2002ll2004l; 
Lueftinger et al. 2003) to reconstruct chemical maps of the sur- 
face of magnetic Ap stars from a series o f hig h resolution spec- 
tral lin e profiles. [Kochukhov et al.l (l2002h and IKochukhov et"a]l 
(12004 demonstrated that even if the surface distribution of some 
chemical species form symmetric patterns closely following the 
magnetic geometry, a multitude of chemical abundance struc- 
tures in the magnetic Ap (and roAp) stars cannot be reduced to a 
system of spots and rings; they also showed that the distribution 
of some elements exhibits clear asymmetry with respect to the 
dipolar magnetic field geometry. 

This suggests that phenomena other than the magnetic field 
interact with microscopic diffusion and play a substancial role in 
determining the geometry of abundance structures (e.g. rotation 
and mass loss). 

2.2.2. Vertical stratification 

Abundances analyses derived from spectroscopic observations 
of magnetic Ap stars have discovered inconsistencies between 
abundances deduced from the weak and strong spectral lines 
and the lines of different ions. These anomalies provide clear 
evidence of prominent vertical abundance stratification of many 
chemical species, such as Ca,Cr, Fe, Na, Mg, Ba , Si, REE el- 
ement s. (lBagnulo etal.!l200ll; IRvabchikova et all l200lL l2002l 



2003t iKochukhovl 120031; IKochukhov et all 12001 IWade eTaL 
2003). 

It appears to be widely accepted that iron (and most iron- 
peak elements) is normal to underabundant in magnetic Ap stars 
upper atmospheres, but numerous pieces of evidence highlight 
the increase in iron abundan ce deeper into the atmosphere as in 
the roAp star HD 213637 (Koc hukhovll2003l) . the roAp star y 
Equ (Rvabchikovaetal. 2002), and the Cr Ap star HD 204411 
(Rvabchikova et al 11200 5). For instance, the stratification model 
of .Rvabchikova et al. f (l2002) for y Equ predicts that Ca, Cr, Fe, 
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Fig. 1. Observations of magnetic Ap stars. Left panel : Hertzsprung-Russell diagram. Solid circles represent 19 roAp stars, open 
circles, 18 noAp stars, and open triangles are 3 roAp star candidates. Right panel : effective temperatures of roAp and noAp s tars and 
roAp s tars candidates. Observations are taken from Kochukhov & Bagnulo (2006) (15 roAp stars + 3 candidates) and .North et al.l 
(Il997h . 



Ba, Si, and Na are concentrated in lower atmospheric layers, but 
are normal to underabundant in the upper layers with a transition 
between -1.5 < log T5000 < -0.5. 

Globally, no obvious differences between the abundance pat- 
tern of Ap stars that pulsate or do not (or are candidate for 
pulsations) have been observed. However, an anomaly is de- 
tected in the REE abundances. In roAp stars, abundances de- 
termined from lines of the third spectrum of various rare earth 
elements (notably Pr III and Nd III) appear to be overabun- 
dant from 1.5 dex and up to 2.5 dex relative to singly ion- 
ized spectra. This anomaly is observed in spectra of all inves- 
tigated ro Ap star s (Weiss et al. 2000; Ryabchikova et al. 2000, 
2001; Gelbmann et al. 2000; C owley et all 120001; IKochukhovl 
2003; .Ryabchikova et al.. .2004), apart from HD 137909 and HD 



1 161 14, which will be discussed later. 

Similar analyses for non-pulsating Ap stars with otherwise 
similar properties did not find evidence of this imbalance o r 
only marginally (IWeiss et all l2000l; iRvabchikova et al.l l2004l) . 
iRyabchikova et alT i 20021) interpreted the roAp star anomaly as 



the result of a stratified Nd and Pr distribution. In the roAp star 
7 Equ, they de tected an acc umulation of these elements above 
logTsooo - -S. lMashonkina et al., (2005.) investigated the forma- 
tion of Nd II and Nd III lines in the atmos phere of y Equ and 
HD 24 712, and confirmed the detection by 'Ryab chikova et al.l 
(I2OO2I) of enhanced Nd abundance layers, although for a differ- 
ent measurement of the transition at logTsooo = -3.5. 

However IRvabchikova et al] (|2002|) did not find a REE 
anomaly in two luminous and evolved roAp stars: HD 137909 
(Teft = 7430K , log(L/Lo) = 1.44) and HD 116114 

i Teff = 7413K, log(L/LH) - 1.32) confirmed to be roAp stars 
Hatzes & MkrtichianI l2004t lElkin et all [2005h . This suggests 
that the true difi'erences between roAp and noAp stars are still 
far from being understood. 



2.2.3. Metallicity of roAp stars 

Significant deviations in the chemical composition of CP stars 
from solar values are believed to be limited to surface layers. 



These anomalies are probably caused by a complex interplay 
between microscopic diffusion (including radiative levitation), 
mass loss, convection, turbulent mixing, and magnetic field. 
These complex effects are not well understood and are difficult 
to simulate in the context of ID stellar evolution computations. 
Therefore, it is impossible to estimate the interior metal content 
from the observed surface abundance patterns. 

According to observations from Gelbmann (1998', Fig. 1), 
Kochukhov (2003, Fig.6) and Ryabchikova et al. (2004, Fig. 5), 
the surface metallicity of the magnetic Ap stars is correlated 
with their effective temperature. The metallicity of stars with ef- 
fective temperatures of between 6000K and lOOOOK increases 
with effective temperature. Iron is underabundant in stars cooler 
than 7000K (by up to one order of magnitude), close to solar 
in stars between 7000K and 8000K, and overabundant in stars 
hotter than 8000K. Since no roAp star hotter than 8700K was 
detected, this could imply that a relationship exists between the 
excitation mechanism and the heavy element distribution. 

3. Theoretical Context 

3.1. Excitation mechanism 

Since the discovery of roAp stars, several mechanisms have been 
proposed and tested extensively to excite their rapid oscilla- 
tions: the direct influence of the Lorentz force ('Dziembowskij 
1984; Dziembowski & Goode 1985), the magnetic overstabil- 
ity (Shibahashi 1983; Cox 1984), the stoc hastic excitation (see 
Houdek e t aL. .199 9), the /^-mec hanism CDolez & Goughlll982l; 
Matthew s) ll98i" Dolez et aQ 119881: lOziembows ki & Goodd 
.1996i:iBalmfor th et al. 2001). For a review of all these processes 
we refer the reader to Balmforth et al. (2001) and Cunha ( 20021) . 
In this paper, we discuss only the most plausible explanation: the 
excitation by the /c-mechanism. 

As a result of their position in the HR diagram, close to the 
5 Scuti stars, it was first suggested that the oscillations of roAp 
stars might be driven by the same mechanism as 6 S c uti modes, 
i Jiat is the /c-mechanism of Hell (e.g. lKurtz|[T990h . iMatthewsl 
(Il988h also suggested that the roAp pulsations could be driven 
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by the A:-mechanism acting in the SilV ionization region, close to 
the Hell ionization region. While no accurate computations were 
carried out to test a possible excitation in the SilV ionization 
zone, several studies tested excitation by the Hell /c-mechanism. 
Numerical computations, h owever, never suppor ted this conjec- 
tur e. As shown byLD ziembo wski & Goodel(ll996l) and confirmed 
bv lBalmforth et al] 02001). only low frequency modes similar to 
those of 6 Scuti stars can be excited in the helium second ioniza- 
tion region. Both sets of authors demonstrated that, in contrast, 
the HI ionization zone can exert significant driving for the high 
order p modes detected in roAp stars. Today the /c-mechanism 
operating in the HI ionization region is considered to be the most 
probable driving mechanism for roAp oscillations. As a result, 
the properties of the excited modes are expected to be strongly 
dependent on the chemical distribution, especially in the external 
layers of the star. 



3.2. Non-standard chemical transport processes 

Several non-standard processes are expected to be involved in 
the definition of the chemical element distribution in magnetic 
A stars. Numerous observations have been found to be consis- 
tent with the hypothesis that microscopic diffusion including ra- 
diative levitation, is a key process in these stars. The inhomoge- 
neous and abnormal distribution of heavy elements observed at 
the surface of the magnetic Ap stars as well as the vertical strat- 
ification detected in their atmospheres (see section 12.2.2b attest 
that element segregation effectively occurs in Ap stars, and that 
competing processes, if any, do not efficiently inhibit its effects. 

Outflows of mass in the form of stellar winds are also 
expected to occur in Ap stars as they are observed in hot- 
ter and cooler main sequence stars. Unfortunately, no accurate 
measurement of mass loss-rates of A stars, only upper limits, 
have bee n available unti l now fr om observations (M < 10 
Mo.yr-i, iLanz & Catalal (Il992l) : iBrown et al.1 (Il990l) . The ob- 
served abundance anomalies in Ap and Am stars may also pro- 
vide an indirect estimation of A-star mass-loss rates: several 
theoretical works demonstrate that mass-loss rates lower than 
10"'^ M0.yr"' (and sometimes as small as 10"'^ Mo.yr"' ) are 
needed to provide quantitative agreement between com puted and 
observed abundances in chemically peculiar A stars (iMichaudl 
[1986: Michaud_eial. 1983; Babd &Michaudl99l). 

Turbule nt transport is hkely to occ ur in stellar radiative in- 
teriors (see ' Zahn & Zinn-JustinI 1 19931 and references therein) 
due to various hydrodynamical instabilities. However, theoret- 
ical works dedicated to the magnetic field/turbulence (or con- 
vection) interaction and observations of abnormal abundances at 
the surface of magnetic Ap stars appear to exclude the presence 
of strong turbulent mixing (see next section) in these stars. 



3.2.1. Magnetic field effects 

The significant magnetic fields detected in magnetic Ap stars are 
expected to influence strongly their oscillation properties. In the 
deep interior, the magnetic field is unlikely to play an important 
role in the dynamics of the oscillations. In the outer layers, it 
is expected, however, to influence the properties of the oscilla- 
tions both directly, by its effects on the mode geometry and fre- 
quencies and indirectly, by its effects on the chemical transport 
processes and the excitation mechanism. 

The magnetic boundary layer influences basic proper- 
ties of the pulsati ons, such as the oscillation frequencies 
and eigenfunctions (iDziembowski & Goodd 119961 : iBigot et al] 



2000i: ICunha & Goughl2000l:ISaio & Gautschvl2004t[Saiol2005l: 
,Cunhall2007h . In particular, it can induce frequency shifts of be- 
tween a few microhertz and a few tens of microhertz. Saio (20051) 
showed that the direct effects of the magnetic field on the oscil- 
lations could stabilize the low radial order pulsations in roAp 
stars, and thus provided an explanation of the absence of these 6 
Scuti type modes in roAp stars otherwise predicted to be excited 
by the stellar evolution standard theory. 

The magnetic field is expected to influence most chemi- 
cal transport processes occurring in A stars. It can deceler- 
ate the microscopic diffusion but this effect is n ot expected 
to be significant in Ap stars ( Theado et alj 12005*). The mag- 
netic field can also prevent stellar winds from flowing close to 
the magnetic equator (see iTheado et alj I2005L and references 
therein). The main effects of the magnetic field, however, are 
related to convection and turbulence. As discussed previously 
(Gough & Tavler 1966; Moss & Tavlor 19691: IShibahashi|[T983i 
fCox 1984; Balmforth et al.. .2001) . a strong magnetic field can 
freeze convective and turbulent motions in Ap stars in a signif- 
icant part of the stars. This conjecture is supported by observa- 
tions (described in Sect. I2.2.2l i which provide clear evidence of 
strong vertical abundance gradients in the atmosphere of mag- 
netic Ap stars, indicating that the magnetic field inhibits convec- 
tive and turbulent mixing in th e outer stellar envelope. 

Observations presented by ' Kochukhov et alj (|2007|) prompt 
a reconsideration of this conclusion. They investigated the vari- 
ations in the REE lines of roAp stars and discovered a change 
in the profile variability pattern with height in the atmospheres 
of all studied stars. They proposed that the line width mod- 
ulation is a consequence of the periodic expansion and com- 
pression of turbulent layers in the upper atmospheres of roAp 
stars. This suggests that some cool Ap stars possess a turbulent 
zone possibly due to convective instability; this would only be 
true, however in the upper atmospheric layers (above log tsooo - 
-3.5), since, in the intermediate atmosphere, significant verti- 
cal gradients of chemical composition preclude substantial mix- 
ing (-3.5 < logTsooo ^ -0.5). This proposition is suppo rted by 
theoretical diffusion models bv' Leblanc & MoriinI (|2004|) . which 
demonstrate that convective mixing is needed to achieve a good 
agreement between stratification profiles inferred from computa- 
tions and observations. These results suggest that the interaction 
between magnetic field and turbulent/convective motions could 
be more complex than previously anticipated. 

3.3. Modeling : status report 

Theorists have attempted to model the complex internal struc- 
ture of magnetic Ap stars to explain, in particular, the properties 
of their oscillations. We now summarize the work completed to 
date. 

Standard theoretical models that adequately describe roAp 
stars predict the excitation of 5-Scuti modes but until now no 
high fre quency modes were found to be unstab le in these mod- 
els (e.g.lDolez & Goughlll982l: IDziembowski & Goodd [19961: 
iBalmforth et alj 1200 lb . As a result, alternative non-standard 
models were proposed to test the influence of various physical 
processes or structures such as a chromosphere, a magnetic field, 
helium settling, stellar winds, or mixing on the excitation mech- 
ani sm of pulsations . 

iGautschv et alj (1 19981) investigated the influence of an ad- 
hoc temperature inversion on the atmosphere of chemically ho- 
mogeneous models. This inversion, which simulates the pres- 
ence of a chromosphere enables roAp pulsations to be excited 
in models of A stars with effective temperatures lower than 
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logTeft - 3.86. This theoretical blue edge for the instability strip 
is far cooler than observed. As emphasized by iGautschv et al.l 
(Il998h . there is also no observational evidence for the presence 
of a chromosphere in roAp stars. 



Fo ll owing work s bvlGough & TavleJ (1 19661) ; iMoss & Tavlod 
(119691) : IShibahashil (Il983l) : Icoxl (|1984|) dedicated to the abil- 
ity of magnetic fields to freeze c onvection, several authors 
jBalmforth et al.ll200ll;ICunhall2002l) tested the stability of mod- 
els with fully radiative envelope. They demonstrated that the 
suppression of convection in the external layers of the stars al- 
lows us to excite roAp pulsations in models with homogeneous 
composition. However, the theoretical instability strip derived by 
citetcunha02 illustrated how these models do not account for the 
position of the roAp stars instability strip: the predicted theoret- 
ical red edge was too hot compared with observations. 



iTheado et al.l (l2005l) presented the most accurate evolution- 
ary models computed to date for roAp stars. They tested the in- 
fluence of helium settling, weak stellar winds, and mild mixing 
on stellar models with convection suppressed in the external lay- 
ers. The stability analyses of some of these models were pre- 
sented in iTheado & Cunha (2006). In models with He settling 
and convection suppressed, the external layers (including the H 
and He ionization regions) become rapidly poor in helium and 
enriched in hydrogen. As a result, excitation by the /c-mechanism 
is expected to be enhanced in the HI ionization region and be 
weaker in the Hell ionization zone. Stability analysis demon- 
strates that diffusion accelerates the damping of 6 Scuti modes 
but the influence of He settling on the excitation of roAp modes 
is insignificant. When a small wind (about lO"''* Mo.yr"') com- 
petes with microscopic diffusion, an accumulation of helium oc- 
curs somewhere below the surface of the star, in the vicinity of 
the helium first ionization region. As expected, the He accumula- 
tion does not affect the excitation of roAp pulsations but favours 
the excitation of 6 Scuti mode pulsations. These models do not 
enable the observed instability strip to be reproduced. In partic- 
ular, they fail to account for the excitation of pulsations in cool 
magnetic Ap stars, i.e. roAp stars observed in the lower right- 
hand-side corner of Fig. [T] 



3.4. Theoretical challenges 

Until now no theoretical study has been able to account for the 
exact position of the roAp instability strip, in particular the ex- 
citation of roAp modes in magnetic Ap stars with effective tem- 
peratures lower than 7400K. On the other hand, the detection of 
both roAp and noAp stars in the upper part of the HR diagram 
and the absence of noAp stars with effective temperatures lower 
than 7550Kboth still present challenges to theorists. Abundance 
determinations in several roAp stars appear to imply that a rela- 
tionship exists between the metallicity and the excitation mech- 
anism. However, no study was dedicated to the influence of the 
heavy element distribution on the pulsations of magnetic Ap 
stars. When examining the most important contributors to the 
opacity in the H ionization zone, we also found, surprisingly, 
that, for surface conditions typical of roAp stars, a main contrib- 
utor is iron, as can be seen in Figs. [T4llA.4l and lB.2l 

The aim of this paper is to study the influence of metals on 
the excitation mechanism of roAp star pulsations and to assess 
whether these elements could help to resolve some open issues 
regarding roAp star properties. 



4. Basic non-adiabatic pulsation equations 

As mentioned previously, roAp stars pulsate in low degree high 
order non-radial p-modes. Since the radial motions of a low de- 
gree high-radial-order-p mode are much larger than the horizon- 
tal motions, the results concerning the stability of radial pulsa- 
tions are applicable to pulsations in roAp stars. We recall the 
basic equations of linear non-adiabatic radial pulsations, since 
they are useful to the forthcoming discussion. The additional 
terms originating in the magnetic field and rotation are disre- 
garded. The time-dependence of the eigenfunction is written as: 
6X{r, f) - 5X(r) exp(/crf). We therefore have 3(cr) < for unsta- 
ble modes. The operator 5 denotes Lagrangian perturbations. 
The perturbed movement equation is given by: 



The continuity equation is: 

So 1 A{r^5r) 
p dr 



(1) 



(2) 



The perturbed equation of energy conservation in the absence of 
nuclear reactions is: 
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The perturbed diffusion equation in a full radiative zone is: 
_ ^ 5r _^ ^ 5r 5k ^ d(ST/T) 

T ~ 7"^ ~ ~ ~ dinr ■ 

The work integral expression is: 



(3) 



(4) 



Sp* 6P\ 
M 



Wim) = - V>{——)Am 
Jo [ P P ) 

(Fs - l)d\—T6s \ dm 
JO [ P 

am 



= -X" 



(5) 



The growth rate of the modes rj {6X{r, t) - 6X{r) exp((i1R(cr) + 
J]) f)) is given by: 



W(M) 



2 %(o-) l^rpdm 



(6) 



The integrant in Eq. (|5]l can be simplified in regions where 
the local wavelength A = 2n/k = Incjcr is far shorter than the 
scale heights of equilibrium quantities. Substituting Eq. ([3]) into 
Eq. (|5]l, retaining only the term d((5r/r)/dlnr in Eq. (HI, and 
neglecting everywhere the eigenfunctions compared with their 
derivatives: 



dW ^ 



|dlnr/dr| 
%{a-)L 



1 51*6^(6711)^ 



!R - — 
a- T 



dr 



c2|dlnr/dr| 



6T 



dr. 



(7) 



Hence the regions of short wavelength oscillations have a 
damping effect, generally refered to as "radiative damping". 
Finally, the perturbed equations of state are: 
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6T 5s „ 6P 



= Pt 



6s 



+ Ti 



dp 



(8) 
(9) 



5. Results 



We computed stellar models to describe A stars using the stel- 
lar evolution code Cles (Code Liegeois d'Evolution Stellaire). 
For a complete descr iption of the code, we refer the reader to 
IScuflaire et alj (|2008|) . The metal mixture used in these compu- 
tations was the solar mixture presented by Asplund et al. (2005|) 
(hereafter AGS05) . The opacity tables are those of OPAL96 
dlglesia s & Rogers! 1 1996 1) completed at low tempera ture with 
tables based upon calculations from [Ferguson et alj (j2005'[[3. 
All t hese tables ar e computed using the composition taken 
from [Asplund et aP (12005). We use the OPALOl equation of 
state jRogers & Navfonov 2002). As outer boundary conditions, 
Kurucz atmospheres CKuruczi il998i) are connected to the in- 
terior at an opt i cal de pth t equal to 1. Follow ing works by 
iBalmforth elal] ([2001) and The ado et all (2005), we computed 
models with and without convection in the envelope. 

5.1. Models with solar metallicity 

5.1.1. Instability strip 

We first consider models with initial mass fractions of hydro- 
gen X=0.71 and solar metallicity according to AGS05 (i.e. 
Z/X=0.0165). These models were computed with a full radia- 
tive envelope. The stability of these models has been computed 
using the non-adiabatic pulsation code MAD (Dupret 2002). We 
restricted our research to the roap frequency domain without in- 
vestigating the possible presence of excited (J-Scuti modes (see 
Sect. I3.2.T] ). The left middle panel of Fig. [TT] shows the evolu- 
tion tracks of our models and the calculated theoretical instabil- 
ity strip of roAp stars. Figure |2] displays, for four masses, the 
frequencies of the unstable roAp type modes along the evolu- 
tionary tracks. Figure |2] shows that, although agreement is good 
down to a luminosity of about log{L/LQ) - 0.9, our models fail 
to reproduce the position of the red edge of the instability strip. 
This is not a surprise since until now no evolutionary computa- 
tion has been able to account for the excitation of roAp modes 
in stars with effective temperatures lower than 7400 K. All com- 
puted models, whether of^ sj;andard homogeneous composition 
jBaTmforth et alj|200ll: ICunhal 120021) or with a st ratified distri- 
bution induced by He settling and stellar winds dTheado et alj 
120051: iTheado & Cunhall2006l) . failed to shift the red Umits to- 
wards cooler temperatures. 

5.1.2. Driving and damping regions 

In the present study, we neglected the coupling of oscillations 
with the magnetic field; for the radiative transfer we use the dif- 
fusion approximation below the photosphere and in the atmo- 
sphere we a dopt the perturbed t reatment of the thermal aspects 
proposed by (Dup ret et al.l2002l) . In this simplified framework, it 
is easier to identify and interpret the main features of the driving 
mechanism in roAp stars. 



' The opacity tables adequate for low temperatures are taken from 
the website http://webs.wichita. edu/phy sics/opacity/ , 
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Fig. 2. Frequencies of the unstable modes predicted by our mod- 
els as a function of log Tetf . The different symbols correspond to 
different masses from 1.6 Mq (cross) to 1.9 Mq (diamond). The 
differences observed in the excited frequency ranges are mainly 
due to the different dynamical times: Tdyn = (i?^/GM)'^^ of the 
models. 



Figure |3]presents the opacities of two 1.6 Mq, [Fe/H]=0.00 
models inside (full curve) and outside (dashed curve) of the the- 
oretical instability strip. In this section we compare in detail the 
results obtained for these two models, represented by stars in 
the left middle panel of Fig.[TT] Figure |4] shows the density dis- 
tribution for these two models. We note the density inversions 
resulting from the large opacity bumps. Figure|5]presents the dif- 
ferential work -dW/d log T for the radial mode p32 (v = 3 mHz 
for the hot model, and v = 1.8 mHz for the cold model). The 
regions where the differential work is positive (negative) have 
a driving (damping) effect on the oscillations. Damping occurs 
between log T = 5 and log T = 4.3 and driving occurs around 
the log r ^ 4.1 opacity bump. 

In Fig. |6] we indicate the contribution of the different terms 
in Eq. |4] The term d(6T /T)/dlnT is given by the dashed line, 
the other contributions to 6L/L are shown by the point-dashed 
line, and the solid line indicates the sum. Two distinct regions 
can be distinguished. From the center to log T ^ 4.5, the 
term d(6T/T)/d In T dominates because the eigenfunctions have 
many nodes there. Substituting this dominating term in Eq. |5] 
shows the well known radiative damping mechanism occurring 
in regions of short-wavelength oscillations. We note that this ra- 
diative damping depends on the equilibrium temperature gradi- 
ent, but not significantly on the fluctuations of the opacity. We 
now consider the region between log T ^ 4.25 and the photo- 
sphere. Because of the strong opacity bump, 6T/T increases dra- 
matically in this region. This can also be seen in Fig. [T] which 
indicates \6T/T\ and \6s/cp\. 

Since the pressure is determined by the movement equation 
ddP/dr = §p(4 -H cr^r/g)(5r/r), it cannot follow the enormous 
temperature variations. For \6P/P\ « \6T/T\, we obtain from 
the perturbed equation of state (Eqs. 8 and 9): Ss/cp ^ 6T/T 
-dpip in the Hi opacity bump (see Fig. |7]i. The heat exchanges 
are therefore extremely large, and significant amounts of en- 
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Fig. 4. Densities of two 1.6 M© models with different Teg. Fig- 6. Contributions to ''R{6L/L] for the radial mode p32 (see 

details in the text), model with Teft = 8090 K. 



ergy can be absorbed from or transferred to the modes. The pre- 
cise result of this driving or damping depends on the respective 
phases of the eigenfunctions. Since 6p/p is more or less in oppo- 
site phase with 6s/cp and 6T/T, it does not help to know if we 
have either driving or damping. It is more helpful to consider the 
phase differences with respect to 6P. From Eq.|5] we can see that 
to achieve a driving of the modes (0° < (p{6P) - (p(6s) < 180°), 
heat must be received {d6Q/dt > 0) when SP > 0. This appears 
to be the case for the radial mode P32 considered here and it is 
therefore excited, but we will see in Sect. 15. 1.4] that the driving 
of a modes strongly depends on the location of its last node. 

The driving mechanism of roAp stars is therefore not identi- 
cal to the weakly nonadiabatic /c-mechanism (as in other clas- 
sical pulsators). In the latter case, the fact that dluK/dlnPl, 
and therefore \6k/k\ increase outwards implies automatically that 
dSL/dm < (Eq. 4), so heat is gained during contraction and 
transformed into mechanical work. In the Hi opacity bump of 
roAp stars, the opacity variations are counterbalanced by the 
temperature gradient variations, and the result of positive or neg- 
ative heat exchanges and mechanical work depend on the modes. 



5.1 .3. What determines the location of the red edge of the 
roAps instability strip ? 

The physical processes determining the instability strip red edge 
differ in our roAp models compared with other /c-driven oscilla- 
tors, such as 6 Scuti stars. In the latter case, the damping of the 
modes origi nates in the time-de pendent coherent coupling with 
convection dPupret et all l2005h . In our purely radiative roAp 
models, we also predict a red edge, and must search for another 
explanation. 

As can be seen in Fig.|3] the Hi opacity bump of cold models 
away from the instability strip is larger than it is inside the insta- 
bility strip. The large temperature gradient in this region implies 
that a density inversion exists that is stronger for the cold model 
(Fig. |4]i. With this in mind, we could have expected more driv- 
ing in the cold model compared with the hot model ( at least in 
a classical K-mechanism vision of the problem) but the converse 
instead occurs. We also note that radiative damping is almost 
identical for the two models (see Fig. [Sj, because the tempera- 
ture gradient and the eigenfunctions are similar in this region for 
different models. The explanation of the red edge comes from a 
closer study of the eigenfunctions and their respective phases in 
the Hi opacity bump. 

As already explained above and shown in Figs. |7] and [8] 
\6slcp\ ^ \6TIT\ « \6plp\ » \6PIP\ in the Hi opacity bump 
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Fig. 7. \6T IT\ and \6slcp\ for the radial mode P32 in two 1.6 M, 
models with different T^g. 
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Fig. 8. \6p/p\ and r^^\6P/P\ for the radial mode P32 in two 1.6 
M0 models with different Teff . 

region. We also note in these figures, that the values obtained for 
\6s/cp\ in this region are close for the hot and cold models. Figure 
[8] illustrates however that the values of \6P/P\ are much smaller 
in the cold model than the hot model. As mentioned before, it is 
more appropriate to consider the work expression in terms of 6P 
instead of dp to be able to observe how much driving occurs in 
this region. 

Since \6P/P\ is far smaller in the cold model, \5(T6s*6P/P)\ 
is smaller and less work is performed (see Eq. |5]l. This smaller 
amount of driving cannot counterbalance the radiative damping 
of deeper layers and all modes are stable. 

The eigenfunctions have a complex non-adiabatic behavior 
in the superficial layers, and the interpretation that we propose 
of the behavior of 6P/P must be considered with caution. A pos- 
sible origin of the differences found for 6P/P between the two 
models could be the densities (Fig.|4|. According to the asymp- 
totic theory, the following behavior of P' /P is expected: 

^(r) oc ^F{jcr/cdr) oc F {Jcr/c dr) , (10) 

where F is a cosine in propagation regions and an exponential 
in evanescent regions, c is the sound speed, and cr the angular 
frequency. The asymptotic theory does not hold of course in the 
highly non-adiabatic superficial layers, but P' /Pis an eigenfunc- 



tion less affected by non-adiabatic effects because of the control 
by the movement equation. At a given local temperature. Fig. 
|4] shows that the cold model has a larger density inversion than 
the hot model, because of its larger temperature gradient. Hence, 
we understand from Eq. ( fTOt that P' /P (and therefore 6P/P) is 
smaller in the cold model. A second complementary explana- 
tion could also be the following. From Eq. (fTOl i, we see that, in 
the asymptotic adiabatic regime, the scale height of the variation 
in the eigenfunctions is given by c/cr, which is acceptable for 
dP/P. However, the scale height of variation in the other thermo- 
dynamic eigenfunctions is -dr/dlnT because they are affected 
by the opacity bump. Since the cold model has a higher opacity 
(Fig. [3]) and therefore larger temperature gradient, the increase in 
6P/P of A In \6P/P\ * (cr/c) Ar * (cr/c) (dr/d In T)AlnT is ex- 
pected to be smaller Hence, 6P/P cannot reach sufficiently high 
values around the Hi opacity bump of the cold model, and high 
order p-modes are stable. 

5.1 .4. Shape and nodes of the eigenfunctions 

We emphasize that the predictions concerning the driving and 
damping mechanisms in roAp stars are extremely sensitive to 
the shape of the eigenfunctions. For these high radial order p- 
modes, nodes can appear in both the driving region and the at- 
mosphere. The latter phenomenon is not purely theoretical, since 
it is observed in the line-profile variations of some roAp stars. As 
an example, we compare in Figs. |9] and [10] the results obtained 
for the same hot model as before, but for two different modes: 
P32 and P37. We can see in Fig. |9] which presents -dW/dlogT, 
that the mode P37 is damped almost everywhere ! To understand 
this result, in the top panel of Fig. [TOl we provide \6p/p\ and 
r^^\6P/P\, and in the bottom panel <p{6p) - (p{6P) (degrees) for 
the modes P32 and P37. We recall (Eq. (|5]l) that regions in which 
0° < 0(5p) - < 180° drive the oscillations, and regions in 

which 180° < (f>{6p)-(f>(6P) < 360° dampen the oscillations. The 
difference between the works shown in Fig.|9]is the presence of 
a node of 6p at log T ^ 3.9 for P37, but not for P32. In both cases, 
(/>{6p) - 4'iSP) - at the surface, because of the small values 
of 6T/T due to the thermal boundary condition. However due to 
the location of the node of dp, 4>(6p) - 4>(6P) enters in the damp- 
ing domain for 4.3 > log T > 3.9, and the mode P37 is therefore 
stable. The location of the nodes generates a type of windowing 
in the driving mechanism of roAp stars. The nodes must be at 
just the right place to allow the driving to occur This window- 
ing explains why the theoretical frequency interval of unstable 
modes is small for any given roAp model (Fig.|2|. 

5.2. Global metallicity influence 
5.2.1. Instability strips 

To test the influence of the global metallicity on the excita- 
tion mechanism of oscillations, we computed grids of models 
for three different metalHcities: [Fe/H]= -0.89, 0.00 and 0.83 
(the initial mass fraction of hydrogen is the same as previously 
X=0.71). Figure nn shows the evolutionary tracks of the com- 
puted models. The left column shows the evolutionary tracks 
of models with convection suppressed and for comparison, the 
right column displays the results obtained for standard models 
(i.e. with convection). The full gray squares represent the mod- 
els for which unstable roAp-type modes are found. 

Figure [12] compares the position of the resulting instability 
strips. The black lines show the position of the instability strips 
derived from models with convection suppressed, and the gray 
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Fig. 11. Evolutionary tracks and theoretical instability strips of model grids with different metallicities. Left column : models with 
convection suppressed in the envelope; right column : standard models (with convection). Upper panel : low metallicity models, 
middle panel : solar metallicity models, lower panel : metal-rich models. The full gray squares represent the models for which roAp- 
type modes are found to be excited. The circles and triangles represent the observations in Fig.[T] For clarity, masses corresponding 
to the evolutionary tracks are given only for a few models. Each model grid is computed with a O.IMq step in mass. The two stars 
on the left middle figure correspond to the I.6M0 models compared in Fig.[3]to[T0l 
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Fig. 9. Comparison of the differential work -dW/d log T for the 
radial modes P32 and P37, 1.6 model with Tgff = 8090 K. 
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Fig. 10. \5plp\ and Y\^\5PIP\ (top panel) and ^{5p) - ^(5F) 
(bottom panel) for the modes ^32 and py], 1.6 model with 
Teff = 8090 K. 



lines the instability strips deduced from standard models. The 
figure illustrates how the suppression of convection increases the 
width of the theoretical instability strip towards cooler temper- 
atures. In contrast, the global metallicity does not influence the 
position of the instability strip significantly. The models with 
solar metallicity, however, provide a closer agreement between 
observational points and theoretical instability strip (see Fig. 
fTTI ). The excitation of modes in stars with effective temperatures 
lower than 7400K (logTeff < 3.87K) is, however, not explained 
by these models. 
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Fig. 12. Superposition of the instability strips deduced from 
models of Fig. [TT] Black lines are the instability strips derived 
from models grids with convection suppressed, gray lines the 
instability strips derived from standard models. 
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Fig. 13. Evolutionary tracks for three models with different 
masses and metallicities showing a common intersection (lo- 
cated by the open black square) in the HR diagram. 



5.2.2. Comparison between calibrated models 

To study in detail the influence of the global metallicity on the 
excitation mechanism of pulsations, we compare three models 
of different masses and metallicities but located at similar posi- 
tion in the HR diagram. Figure[T3]shows the evolutionary tracks 



of the three computed sequences. In the following, we compare 
the structures of the models located at the intersection of the 
evolutionary tracks. The characteristics of the three models are 
presented in Table [T] 

Figure[T4]compares the opacity distribution in the three mod- 
els. The global metallicity strongly affects the opacity profiles. 
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Table 1. Global characteristics of three models with different 
global metallicities and the corresponding frequency range of 
excited modes. 



M(Mo) 


[Fe/H] 


Teff (K) 


log(L/Lo) 


Excitation range 


1.3994 


-0.89 


7715 


1.3919 


0.666-0.698 mHz 


1.8358 


0.00 


7690 


1.3925 


0.902-1.009 mHz 


2.8000 


0.83 


7714 


1.3918 


1.526-1.569 mHz 



150 



1.3994Mq [Fe/H] = -0.89 
1.8358Mq [Fe/H] = 0.00 
2.8000Mq [Fe/H] = 0.83 



I ; 




Fig. 14. Opacity profiles in models with different masses and 
metallicities but located at the same point in the HR diagram. 



In particular, it affects the opacity close to logT ^4.1, i.e. in the 
H ionization zone. In the context of a "normal" /c-mechanism 
driving, we should expect the metallicity to affect the excitation 
of the modes. 

However, this does not appear to be the case, since the three 
models have excited high order p-modes. To understand more 
clearly what is happening, we provide in Fig.[T5]the differential 
work obtained for these three models. The modes considered in 
this figure are: 

- for the [Fe/H]=-0.89 model: p2i (v = 0.698 mHz, excited), 

- for the [Fe/H] =0.00 model: P25 (v = 0.938 mHz, excited) 
and P21 (v = 0.792 mHz, stable), 

- for the [Fe/H]=0.83 model: P34 (v = 1.569 mHz, excited) 
and P21 (v = 0.979 mHz, stable). 

We can see in Fig. [15] that the driving is similar for the 
three excited modes, although they have significantly different 
radial orders. On the other hand, no significant driving occurs 
for the P21 mode of the [Fe/ H] =0.0 and [Fe/H] =0.83 mod- 
els. As we explained in Sect. 15.1.41 a key point for the driv- 
ing of roAp modes is the last node location of the eigenfunc- 
tion \6P/P\. In Fig. [16] we show \6P/P\ for the same modes and 
models as before. We can see that the last node of \SP/P\ is lo- 
cated at the same temperature for all three excited modes. To 
have efficient driving, it must be located at log Ti ^ 4.35, just 
between the Hi and Hen opacity bumps. The location of this 
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Fig. 15. Comparison of the differential work -dW/dlogT for 
three models with different global metallicities ([Fe/H]=-0.89, 
[Fe/H]=0.00 and [Fe/H]=0.83) at the same location in the HR 
diagram. 



last node corresponds to a specific value of the acoustic depth 
given by = J cr/c dr - J^ ' cr / c (dr / dT) dT , where Ti is the 
temperature at the last node, and T2 is the surface temperature. 
Conversely, from the relation: 



o- = T, / c-\dr/dT)dT 



(11) 



we can infer the frequency of the mode whose last node is lo- 
cated at a fixed temperature (Ti) and a fixed acoutic depth (tc). 

This simple relation allows us to understand how changes 
in the metallicity of the models change the range of unstable 
modes. Increasing Z increases the opacity significantly as shown 
by Fig. [14] and therefore the temperature gradient. Hence, cr 
must increase, according to Eq. (fTTT i. In other words, as dT/dr 
increases, the size of the opacity bump region decreases, so that 
we have to decrease the wavelength by increasing the frequency 
to ensure that the last node remains at the same temperature. In 
Fig. [16] the location of the last node of the p2i modes confirms 
that this reasoning is correct. If the frequency of the mode is not 
increased sufficiently, the size of the opacity bump region de- 
creases more rapidly than the wavelength so that the last node is 
located at higher temperature (at log T ^ 4.4, dotted curve for 
the [Fe/H] =0.00 model, and at log T ^ 4.5, dot-dashed curve for 
the [Fe/H]=0.83 model), which inhibits the driving. 

Since the location of the last node is the critical point for 
the driving, and since it is always possible to find a mode with 
this location, we can understand why the location of the theoret- 
ical instability strip is poorly sensitive to changes in the structure 
model, even if they strongly affect the opacity bump in the driv- 
ing region. 

6. Conclusions 

Observations appear to imply that a relationship exists between 
the excitation mechanism of pulsations in roAp stars and their 
heavy element distribution (cf Sect. \2.2.2i . On the other hand, 
the region that plays a major role in driving the roAp high 
order p-modes corresponds to the last opacity bump around 
logTeff - 4.1. It has been assumed that this opacity bump was 
due mainly to the partial ionization of hydrogen. However, our 
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Fig. 16. Comparison of \6P/P\ for three models with dif- 
ferent global metallicities ([Fe/H]=-0.89, [Fe/H]=0.00 and 
[Fe/H]=0.83) at the same location in the HR diagram. 



study has demonstrated that the contribution of metals (typically 
iron) is significant. We have therefore determined how the theo- 
retical instability strips of roAp stars is affected by a change in 
metallicity. Models with different global metallicity (Sect. 15.21 ). 
local enhancement of metals (Appendix |A] available in elec- 
tronic form only), or iron (Appendix |B] available in electronic 
form only), were considered. The local enhancements are ex- 
pected to be produced by the competing effects of gravitational 
settling and radiative levitation. Although the last opacity bump 
of these models is affected significantly by these changes, it ap- 
pears that the theoretical instability strips are only poorly af- 
fected. Our models reproduce well the blue edge of the roAp 
instability strip, although the red edge deduced from the models 
is always too hot, so that our instability strips do not account for 
the cold roAp stars. 

We find that the driving mechanism of roAp high order 
p-modes differs qualitatively from the weakly nonadiabatic k- 
mechanism acting in other classical pulsators. The shape of the 
eigenfunctions and, in particular, the location of their last node 
is critical. To allow efficient driving to occur, the last node of 
6P/P must be located between the two last opacity bumps at 
logT ^ 4.3. Changes in the structure of the superficial layers 
affect the driving less than the eigenfunctions. It follows that 
changing the opacity affects the driving differently than the clas- 
sical /f-mechanism. In some cases, what we could refer to as 
an "inverse ^--mechanism" can occur: increasing the opacity and 
thus the temperature gradient leads to a density inversion close 
to the photosphere; the larger temperature gradient and density 
lead to smaller SP/P and therefore a smaller amount of driving, 
which produces exactly the opposite effects of the classical k- 
mechanism. This explains why the theoretical instability strips 
are less affected by the changes in our models than expected on 
the basis of the classical A:-mechanism. 

The weak effect of a drastic increase in opacity in the driv- 
ing region is a surprising result, which led us to investigate the 
details in the A:-mechanism at work in such stars. In the future, 
we propose to consider additional changes to the models and 
observe how they affect the driving of the roAp high order p- 
modes. In particular, we intend to consider in more detail local 
changes in the He abundance of different regions (e.g. Hell and 
Hel ionization regions). 
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Appendix A: Local metallicity enhancement 

As described in Sect. 12.2.21 there is observational evidence for 
vertical stratification in cool magnetic Ap stars, and for ac- 
cumulation of some heavy elements in the lower atmosphere. 
Moreover, the coolest magnetic Ap stars have underabundances 
in some iron peak elements, which could also be an indication of 
accumulation of these elements somewhere below the surface. 

Diffusion calculations in model atmospheres by 

iLebl anc & Monin (2004) suggested that, in cool Ap stars 
a high concentration of iron-peak elements occurs in low 
atmospheric layers and that it is plausible to expect a sharp 
decrease in Fe abundance above an optical depth log T5000 - - 1 . 

The present version of Cles does not allow us to compute 
models for A stars that include consistent computations of ra- 
diative acceleration effects. We therefore model the effects of 
microscopic diffusion in a parametric way. In A stars, micro- 
scopic diffusion is expected to lead to accumulation of metals in 
some regions and to metal depletion in others. To locate these 
regions, we computed the radiative accelerations on Fe and Ni 
in several ZAMS models. The radiative accelerations were com- 
puted using the OPCD_2.1 package taken from the websit^^of 
the Opacity Project (see Seaton 2005, and references therein) . 
The OPCD_2.1 package contains data, codes and instructions for 
computing Rosseland mean opacities and radiative accelerations 
for any given stellar structure. Using this package, we computed 
the radiative acceleration on iron and nickel in ZAMS models 
with masses between 1.4 and 2.OM0. The initial hydrogen mass 
fraction of these models was X=0.71 and the metallicity was of 
the solar value. The upper panels of Fig. lA. II compare the ra- 
diative accelerations on iron and nickel to the local gravity in 
1.4, 1.6, 1.8, and 2.0 Mq ZAMS models. The lower panels show 
the derivatives of the radiative accelerations. Iron and nickel are 
expected to accumulate in regions where the derivatives of the 
radiative accelerations are positive: we therefore expect iron to 
accumulate close to log T - 4.05, 4.4, 4.6, and 5.2, and nickel to 
accumulate just below the surface (down to log T = 4.1), close to 
log T = 4.35, 4.6 and 5.4, at least during the early main-sequence 
phase. 

To study the influence of diffusion-induced iron/nickel accu- 
mulations on the excitation mechanism of magnetic Ap stars, we 
adopted an approach of applying the Fe and Ni enhancements 
to the entire mixture. We computed models in which paramet- 
ric metal profiles, centered on chosen temperatures were intro- 
duced to simulate, qualitatively, iron peak element accumula- 
tions. Since the relevant regions for the driving of pulsations are 
the external layers where hydrogen and helium experience their 
ionization, we chose to study the influence of metal accumula- 
tion down to logT = 4.6. 

Our parametrization of the metal profiles was inspired by Eq. 
4 of Balmforth et al. (2001). To center the metal accumulation 
profile on the HI ionization region, we used the following for- 
mula : 

Z = Z, + A[4xi2(l - xn)V + xniZi - Z,) 

where the profile centered on the Hel ionization region obeys the 
law : 

Z = Z, + A[4x2i(l - X2i)]^ + (1 - X2i){Zi - Zs) 

and the profile centered on the Hell ionization region results 
from: 

Z = Z, -H A[4X23(1 - X23)]-^ + X23(Zi - Z,) 
^ http://vizier.u-strasbg.fr/topbase/op.html 
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Fig. A.2. Metal accumulation profiles (mass fraction) introduced 
in 1.7Mq models. The model labelled pZl present a metal accu- 
mulation centered on the HI ionization region. In the pZ2 and 
pZ5 models the accumulation is centered respectively on the 
Hel and Hell accumulation region. The pZ3 and pZ4 models 
present an accumulation centered respectively on logT = 4.35 
and log T = 4.40. 



where x,y is the ionization fraction of element i (i being its nu- 
clear charge) in the (j-l)-th state of ionization. As an example, 
X23 is the ionization fraction of doubly ionized helium. The pa- 
rameters A and Z^ are chosen arbitrarily to equal 0.075 and 
0.65x0.01 17 respectively. To obtain profiles centered on a given 
temperature (different from those of HI, Hel and Hell ionization 
regions), we use one of the aforementioned parametrization but 
we shift the profile towards the chosen temperature. Figure |A!2] 
displays the profiles introduced into the various computed grids 
of models. The results are shown only for the 1.7Mq models; the 
profiles introduced into the other masses originate in the same 
parametrizations and are similar. 

Figure |A3]presents the theoretical instability strips obtained 
for grids of models including the metal parametric accumulation 
profiles of Fig. IA.2I 

We see that here also, the theoretical instability strips are 
not affected significantly by changes in the models. To under- 
stand more closely what is happening here, we consider four 
models at the same location in the HR diagram: M = 1.8 M©, 
Teif = 7600 K, log(L/Lodot) = 1-209, with diff^erent Z profiles: 
constant, pZl, pZ2, and pZ5, according to Fig. IA.2I These four 
models are located close to the red edge of the instability strip. 
We present in Fig. |A.4| the opacities found for these four models. 
The opacity bumps are higher in the region where Z is increased. 
Figure |A3] shows the differential work found for the mode p2y in 
these four models. Despite the differences in the opacity and for 
the same reasons discussed in Sect. 15.1.31 and 15.2.21 no signifi- 
cant difference in the driving mechanism is found between these 
four models, which explains why the location of the theoretical 
instability strip is not significantly affected by local changes in 
Z. 
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Fig. A.l. Upper panels : radiative accelerations on iron and nickel versus temperature inside four ZAMS models with different 
masses. The nearly horizontal curves represents the local gravity. Lower panels : derivatives of the radiative accelerations on iron 
and nickel. These models are computed with X=0.71 and [Fe/H]=0.00. 



Appendix B: Iron local enhancement 

In this section, we investigate the effects of heavy element accu- 
mulation with a different approach. We study models including 
an iron local enhancement but a constant metal mass fraction 
in the radiat ive interior. These co mputations are motivated by 
results from iRichard et alj (12001 1) , which demonstrated that an 
iron accumulation due to radiative levitation may occur in the 
external regions of A stars while keeping the metal mass frac- 
tion constant. 

We introduce in our models a parametric iron profile to sim- 
ulate qualitatively a diffusion-induced iron accumulation. We 
center the iron accumulation profile on the HI ionization region 
where most of the driving is supposed to take place, and where 
iron is expected to accumulate. The metal mass fraction remains 
constant in the radiative interior during the entire main-sequence 



evolution, but the composition of the metal mixture varies with 
depth: the iron abundance is increased in the external layers, 
while the other elements are decreased proportionally, keeping 
the metal mass fraction Z constant. The initial hydrogen mass 
fraction is the same as previously (X=0.71), the metallicity is 
still the solar one (Z/X=0.0165, AGS05) and the models have 
fully radiative envelope. 

In these computations, the opacity is computed by interpola- 
tion between five tables built with different relative abundances 
of iron in the metal mixture (FexO, Fexl, Fex2, Fex5, FexlO), 
the other elements being decreased proportionally). These ta- 
bles were computed from the OPAL website; they were com- 
pleted at low temperature with tables computed specifically by 
J. Ferg uson and based upon calculations from iFerguson et al.l 
(l2005l) . 
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Fig. A.4. Opacity in models with different profiles of Z (constant. Fig. B.l. Iron accumulation profiles introduced in a 1.7Mq 

pZl, pZ2 and pZ5 according to Fig. lA.2l ) but located at the same model, 
point in the HR diagram. 
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Fig.A.5. Comparison of the differential work -dW/dlogT in 
models with different profiles of Z (constant, pZl, pZ2 and pZ5 
according to Fig. |A.2t but located at the same point in the HR 
diagram. 



As an example, Fig. lB.ll displavs the iron accumulation pro- 
file introduced in the I.VMq model. The same profile is used for 
all the computed models with masses between 1 .3 to 2.2M0. The 
resulting opacity profile (referred as "pFe") is shown in Fig. IB. 21 
The "solar" and "pZl" labelled curves represent the opacity pro- 
files in two other 1 .7Mo models located at the same place in the 
HR diagram (Teff=7792K, log(L/Lo)=1.08) near the red edge of 
the theoretical instability strip: the "solar" model includes the 
solar metal mixture (AGS05), the pZl-model includes the metal 
accumulation profile as shown in Fig. lA.2l 

Figures IA.2I IB. 11 and IB. 21 show that increasing (in an ade- 
quate way) the iron mass fraction while keeping Z constant may 
produce effects on the opacity similar to a local metallicity in- 
crease. This shows that the opacity in the HI ionization region 




Fig.B.2. Opacity in a main sequence 1.7Mo model (labelled 
"pFe") including the iron accumulation profile displayed in Fig. 
IB. II For comparison the opacity of two others I.7M0 models lo- 
cated at the same point in the HR diagram are also shown: the 
model labelled "solar" includes the AGS05 solar metal mixture, 
the pZl model is the one displayed on Fig. lA.2l 



is indeed sensitive to the iron abundance, which could therefore 
affect the excitation mechanism of Ap stars. 

The right lower HR diagram in Fig. |A3] displays the theo- 
retical instability strip deduced from this new grid of models. 
This Fig. together with the middle left panel of Fig. [TT] show 
that an iron local enhancement around logT=4.05 leads to a the- 
oretical roAp star instability strip similar to that deduced from 
chemically homogenous models of solar composition. As a con- 
sequence, diffusion-induced iron or metal accumulations in the 
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external layers of magnetic A stars do not account for the posi- 
tion of the observed instability strip of roAp stars. 
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Fig. A.3. Figure pZl to pZ5 : theoretical instability strips deduced from grids of models including the metal accumulation profiles 
shown in Fig. IA.2I Figure pFe : theoretical instability strip deduced from models including the iron accumulation profile presented 



